Solving the Coincidence Problem: Tracking Oscillating Energy 
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Recent cosmological observations strongly suggest that the universe is dominated by an unknown 
form of energy with negative pressure. Why is this dark energy density of order the critical density 
today? We propose that the dark energy has periodically dominated in the past so that its pre- 
ponderance today is natural. We illustrate this paradigm with a model potential and show that its 
predictions are consistent with all observations. 



Introduction. A variety of evidence accumulated over the 
last several years points to the existence of an unknown, 
unclumped form of energy in the Universe. First was 
an apparent concordance Q of different measurements: 
the age of the Universe; the Hubble constant; the baryon 
fraction in clusters; and the shape of the galactic power 
spectrum. Second came the stunning observations ||] 
of tens of distant Type la Supernovae, which found a 
distance-redshift relation in accord with a cosmological 
constant, but in strong disagreement with a matter dom- 
inated Universe. Finally, this past year has seen analy- 
ses H of the experiments measuring anisotropics in the 
CMB. Taken together, the CMB experiments plot out a 
rough shape for the power spectrum, one that is in accord 
with a flat Universe, but in disagreement with an open 
Universe. If we believe the estimates of matter density 
coming from observations of clusters Q], the only way 
to get a fiat Universe, and hence account for the CMB 
measurements, is to have an unclumped form of energy 
density pervading the Universe. 

Perhaps the simplest explanation of these data is that 
the unclumped form of energy density corresponds to a 
positive cosmological constant j^]. A non-zero but tiny 
constant vacuum energy density (cosmological constant) 
could conceivably be explained by some unknown string 
theory symmetry (that sets the vacuum energy density 
to zero) being broken by a small amount. However, to 
explain in this way a constant vacuum energy density 
of 2 x 10~ 59 TeV 4 , which is not only small but is also 
just the right value that it is just beginning to dominate 
the energy density of the Universe now, would require 
an unbelievable coincidence. A different possibility is to 
give up the dream of finding a mechanism which would 
set the vacuum energy density to exactly zero and resort 
to believing that anthropic considerations select amongst 
> igioo s t rm g vacua to find one with a vacuum energy 
density sufficiently fine-tuned for life. Although this an- 
thropic selection mechanism is logically consistent and 
even predicts a small but observable cosmological con- 
stant, one might think that nature would have found a 
more efficient mechanism to obtain a sufficiently small 
cosmological constant than such extreme brute force ap- 
plication of anthropic selection. 

An alternative is to assume that the true vacuum en- 



ergy density is zero, and to work with the idea that the 
unknown, unclumped energy is due to a scalar field <fi 
which has not yet reached its ground state. This idea, 
which is called dynamical lambda or quintessence, has 
received much attention || over the last several years. 
However, two problems still remain. First, the field's 
mass has to be extremely small, less than or of order 
the Hubble constant today ~ 10 _33 eV, to ensure that 
it is still rolling to its vacuum configuration. This is 
in general difficult because scalar fields tend to acquire 
masses greater than or of order the scale of supersym- 
metry breaking suppressed by at most the Planck scale: 
to > F/mpi > TeV 2 /m Pl ~ 1(T 3 eV. Although difficult, 
this could be achieved using pseudo-Nambu-Goldstone 
bosons 0. Another more speculative way to achieve 
this would be to use the hypothetical symmetry (per- 
haps some sort of hidden supersymmetry) that ensures 
that the true vacuum energy density is zero to also pro- 
tect the flat directions in scalar field space that would 
correspond to the very light scalar fields necessary for 
quintessence. The second, and perhaps even more serious 
problem is that almost all of these models require that 
we live in a special epoch today, when the quintessence is 
just starting to dominate the energy density of the Uni- 
verse, and furthermore this specialness cannot even be 
justified by use of anthropic arguments. 

In recent years a lot of progress has been made in un- 
derstanding the behavior of quintessence fields. A broad 
class of solutions, called tracker solutions |10|, has been 
discovered in which the final value of the quintessence 
energy density is insensitive to the initial conditions. For 
example, potentials like V — Vo(f>~ n or V = Vq exp(l/<^) 
can, for suitable choices of Vq, catch up with the crit- 
ical density late in the evolution of the Universe for a 
wide range of initial conditions and thus provide a natu- 
ral setting for explaining the current acceleration of the 
Universe. However, the suitable choice of Vb must be of 
the order of the critical energy density today, i.e., we are 
back to the problem of living at a special epoch today 
and not even being able to use anthropic arguments to 
justify this specialness. 

In a subset of these tracking models, which we call 
the exact tracker solutions [||J|], the scalar field energy 
density is always related to the ambient energy density in 
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the Universe: if the dominant component in the Universe 
is radiation, then the tracking field's energy density also 
falls off as a -4 , where a is the scale factor of the Universe. 
If the dominant component is matter, then the field's en- 
ergy density scales as et~ 3 . This behavior arises from 
an exponential potential for (f> (regardless of the value 
of Vo). Since the energy density in this field is always 
comparable to the background density, we are not living 
at a special epoch: any observer in the distant past or 
future would also see the tracking field's energy density. 
However, these tracking solutions run into two problems. 
First, if their energy density today truly is dominant, 
then it should also have been dominant at the time of 
Big Bang Nucleosynthesis (BBN). Constraints from ob- 
servations of light element abundances preclude such an 
additional form of energy density at early times. Sec- 
ond, tracking models have the wrong equation of state at 
present since the tracking field behaves like matter, with 
zero pressure, instead of having the necessary negative 
pressure to accelerate the Universe. 

In this letter we ask the question, what if the Universe 
has been accelerating periodically in the past? Then the 
fact that the Universe is accelerating today would not be 
surprising. It would merely reflect that the period is such 
that the Universe is accelerating today. Of course, if it 
turned out that to achieve a presently accelerating Uni- 
verse the period had to be excessively fine-tuned, then 
this scenario would not be worth considering. However, 
note that the assumption that there is nothing special 
about the present time itself argues for the robustness 
of such a scenario. If the Universe does accelerate pe- 
riodically, then there is no reason why it should not be 
accelerating today. If the Universe does accelerate pe- 
riodically, then it is, in fact, reasonable to expect it to 
accelerate today. 

To judge the merits of this scenario in a concrete man- 
ner, we adopt an ad-hoc potential. Though worked out 
for this specific potential, the predictions outlined here 
are the generic predictions of a periodically accelerating 
Universe. The model we adopt for study is a modifica- 
tion of the exponential potential (which leads to the exact 
tracker solution). The modification to the potential is a 
sinusoidal modulation, which induces the tracker field to 
oscillate about the ambient energy density. We show that 
such a potential can satisfy the the BBN constraints, can 
produce the right equation of state today and leads to 
testable features in the CMB and matter power spectra. 
We call this type of energy Tracking, Oscillating Energy, 
or TOE. 

The potential and the field evolution. Consider a scalar 
field 4> with potential V((f>) — Vo exp(— A</)v8ttG). It is 
well-known || that such a potential leads to an attractor 
solution with f2^ = p^Kp^f, + p ) = n/X 2 where p Q is the 
energy density in the other component of the Universe, 
which is assumed to scale as a~ n . Thus, no matter what 
the initial conditions are for <fi, it always evolves so that 



it tracks the rest of the density in the Universe. 
Now consider the potential 



V(</>) = Vo exp 
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1 + Asm(u4>VSTrG) . (1) 



This potential serves to modulate the tracking behavior. 
Figure |l| shows the resultant evolution of cf> and its en- 
ergy density for a particular set of the parameters A, v. 
(The normalization Vq can be set to G~ 2 by shifting the 
initial value of </>.) Also shown is the tracking solution 
for this particular value of A without the modulation. As 
expected, the sinusoidal term in the potential leads to 
oscillations about this tracking behavior. One can ob- 
tain analytic solutions for the dynamics of the potential 
in Eq. (1) during radiation (n — 4) or matter (n = 3) 
domination in the limit that A is small by perturbing 
about the corresponding exact tracker model which has 



inG = j In a. The sine in Eq. (1) provides a periodic 
forcing term with period In a = , while the natural pe- 
riod |s) of the damped oscillations about the exact tracker 
solution is In a = 8tt\/ — n)[3(37i — 2) A 2 — 8n 2 ] with 
decay e-life In a = 4/(6 — n). Although the above re- 
sults are strictly valid only for small A, they account re- 
markably well for the behaviour shown in Figure 1. The 
forced period corresponds to the longer period of 5.4 units 
(n = 4) and somehwere between 5.4 units and 7.1 units 
(n = 3), while the natural period corresponds to the 
shorter period of 1.6 units (n = 4,3) of the damped os- 
cillations which are presumably excited by the non- linear 
effects that appear when A is not small. 



1 = 4.0, A = 0.98, v = 0.51 
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FIG. 1. The fraction of the critical density in <f> for the 
potential in Eq. ([!]). The dotted line shows the corresponding 
tracking solution (A = 0) . The upper set of curves shows the 
evolution in (j> for the TOE and the tracking models. 

The energy density due to <p is relatively small at the 
time of BBN and relatively large today for the parame- 
ter set in Figure |l|. It is, of course, clear that in order 
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to get the right behavior at BBN and today, one has to 
pick the "correct" parameter sets. This involves a bit of 
fine-tuning which, as we argue below, is quite reasonable 
and natural. If one thinks of the parameter set as being 
randomly selected, then there is a finite probability that 
the Universe will be accelerating today and that the en- 
ergy density of <fr will be sub-dominant at BBN. What 
is this probability? If one selects A, v and A randomly, 
the chance of getting a Universe like ours is of the or- 
der of 1 in a 100. The exact number (for this potential) 
depends on how stringently we define "a Universe like 
ours". For example the tight constraints 0.4 < Q<f, < 0.8, 

< —0.5, and (p<p / Po)bbn < 0.1 give a probability of 1 
in 450, while the relaxed constraints 0.1 < Q$ < 0.9 and 

< —0.25 and (p<p/ Po)bbn < 0.2 give a probability of 1 
in 26. It is also very important to note that whatever the 
extent of fine-tuning, all of it is in dimensionless numbers. 
There are no energy scales in this scenario which are to 
be set by the present expansion rate of the Universe. 
Power Spectra. To compare with CMB and large scale 
structure observations, we compute the power spectra 
of the perturbations in a TOE model. Perturbations 
evolve differently in the presence of the scalar field en- 
ergy density. For example, perturbations typically grow 
only when the Universe is matter dominated. Therefore, 
we expect a non-zero £l<p to lead directly to power sup- 
pression on the scales inside the horizon, with increased 
suppression for larger O^. 
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FIG. 2. The angular photon power spectrum from the TOE 
model of Figure |y. Also shown is a cosmological constant 
model with all other parameters equal. 

The prediction for the CMB angular power spectrum 
is plotted in Figure ||. The primeval power spectrum 
is scale-invariant with adiabatic initial conditions. Also 
plotted for comparison is a model (ACDM) with cosmo- 
logical constant Qa = Q<j> today and the rest of the cos- 
mological parameters also being the same. In further dis- 
cussions we will contrast the results from the TOE model 



against this ACDM model. A noteworthy feature in Fig- 
ure H is the increase in the heights of the first two peaks 
compared to that of the ACDM model. This stems from 
the fact that the gravitational potential decays more in 
the presence of the additional quintessence energy den- 
sity. The decay of the potential at and after recombi- 
nation (the so-called Integrated Sachs- Wolfe , or ISW, 
effect) leads [[l2| to enhanced power on scales I < 600, 
after which the potential becomes irrelevant. Note that 
the increase in the amplitude of both the first and second 
peak cannot be mimicked by adding more baryons, which 
raise the odd peaks but lower the even ones. 

On smaller scales (I > 600), the TOE model has 
smaller anisotropies. Here there are two competing ef- 
fects. First, the difference between the TOE and the 
ACDM models (around recombination when A is insignif- 
icant) is the presence of the extra quintessence energy 
density, which leads to the expansion rate in the two 
models being related as- 

ffroE(a) = ffACDM(a) x (1 - fy(a)) -1/a . (2) 

Eq. H implies that all the relevant scales at recombina- 
tion (which occurs at a r ~ 10 -3 ) are smaller in the TOE 
model by a factor of about */l — n^(a r ). In particu- 
lar, the damping scale is smaller, which increases in the 
power on small scales for the TOE model relative to the 
ACDM model. The second effect is the large scale nor- 
malization of the two models jljj , and this second effect 
more than compensates for the first. COBE normaliza- 
tion is sensitive to scales around I = 10 for which the dif- 
ferences in the two models with regard to the late-ISW 
effect is important. In particular, since A domination 
occurs very late, the ISW contribution around £ = 10 
is much larger in the TOE model. This in turn implies 
that the normalization of the primeval power spectrum 
is smaller, a fact noticeable in the smaller amplitude of 
the photon power spectrum for the TOE model at small 
scales (and also the matter power spectrum, as we will 
soon see) . One last effect that is worth pointing out con- 
cerns the difference in the peak positions in the two mod- 
els (though unlike the peak amplitudes, it is probably not 
easily discerned). In particular, the TOE model has the 
acoustic features in its angular power spectrum shifted to 
smaller scales. This directly traces to the decrease in the 
angular diameter distance to the last scattering surface, 
for the TOE model. Of course, there is also the compet- 
ing effect of the decrease in the size of the sound horizon 
at last scattering for the TOE model, which minimizes 
the effect. 
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FIG. 3. The matter power spectrum from the TOE model 
of Figure [l]. Also shown is a cosmological constant model with 
all other parameters equal. Power is significantly smaller in 
the TOE model. 

The prediction for the matter power spectrum is plot- 
ted in Figure |^. The difference in power at the largest 
scales is due to COBE normalization and the difference 
in the super-horizon growth factor (which is sensitive to 
the equation of state of the cosmic fluid) for the pertur- 
bation. As one moves to smaller scales, which entered 
the horizon well before the present, the differences in the 
evolution of the matter perturbation become more pro- 
nounced. The presence of the extra quintessence energy 
stunts the growth of perturbation once a mode enters 
the horizon. So, the earlier the mode enters the hori- 
zon, the larger the growth suppression relative to the 
ACDM model. In other words, smaller modes are mono- 
tonically more suppressed (something that may not be 
noticeable in the log plot) compared to the same modes 
in ACDM model. It might also be surprising that the 
4> domination around a — 10~ 6 does not cause a more 
appreciable feature (i.e., suppression) in the power spec- 
trum. The reason is that the smallest scales in Figure || 
have just entered the horizon at the time of <$> domination 
(a - 10~ 6 ). 

The normalization on the small scales is generally 
quoted in terms of as, the rms mass fluctuation within 
a8/i _1 Mpc sphere. For the parameters in Figure [I], the 
TOE model has as — 0.4. This is several sigma smaller 
than the preferred value (see e.g. [jll]]) of ~ 0.8, but could 
be rectified by a small blue-shift in the primordial spec- 
trum @. 

Conclusions. We have constructed a model wherein the 
energy density tracks the dominant component in the 
Universe; satisfies the BBN constraints; and has the 
proper equation of state today. Further, this model 
makes definite predictions for large scale structure and 
for the CMB. 

Perhaps the greatest drawback of this class of models is 



the arbitrariness of the potential. In particular we know 
of no theory which predicts a potential of the form given 
in Eq. ([!]). Nonetheless, we feel that the testable predic- 
tions of the model and the aesthetic quality it preserves 
that we do not live in a special epoch are of sufficient 
interest to warrant further study. 
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